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ABSTBilCT 


The intonation of a speech is an indication 
of the type of statement of the sentence eg, an 
interrogative, a declarative etc. Many important 
inference can he deduced hy knowing the intonation 
pattern in a continuous spee'ch recognition system. 
Intonation are of two types: word intonation and 
sentence intonation. A method has been developed 
to find out the sentence intonation here. 

Vibrations of vocal cords result in the fundar- 
ments,l frequency (f^) of speech. The variation of f 
with time is the intonation contour. 

The variation of f^ has been evaluated with 
respect to time by digital computation. The e^nalog 
speech signal after conversion to digital form is 
processed throu,,h a bank of band pvoss filters simulated 
in the digital computer. The output of this bank is 
stored in a C-matrix. The contents of C-matrix indicate 
the intonation pattern of the utterance. 
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CHAPTER 1 


IRTROPUGTIOH 

1.1 Motivation 

With the vast development in the digital computers, 
the communications with machines have already made the 
transition from science fiction to a practical reality. 
Modern computers now can handle problems considered to 
be non-trivial and complex by h\iman beings. This has 
opened the path for more versatile, ’natural man-computer 
communication techniques. Thus recent years have seen the 
progress of time sharing systems for rapid man-machine 
interaction! graphical display of pictures, graphs, 
languages on cathode ray tubes; and is discussed herein 
an approach for man-machine commimication by speech(voice) . 
Among the various channels available for communication with 
machine, the speech has an edge over the others [ 23 ] • 

A few advantages of speech channel are; it 
provides most effortless encoding of all output channels, 
has higher data rate compared to other output channels, 
prefarable for spontaneous output, does not tie up hands, 
eyes, feet or ears; can be used while in motion, requires 
inexpensive and readily available terminal equipment. 

Since the future man-machine communication would 
be desirable through voice, dus to the above advantages 
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and since a comp-uter is a most versatile machine we are 
interested here to findout how digital computers can he 
used in speech recognition systems. 

1.2 Role of Comnuter in Speech Recognition 

The speech recognition hy computers can he viewed 
from two angles. 

(a) The contents of the message 

(h) The linguistic aspects of the message. 

In the former, the words used in the speech are analysed 
to gather certain information for the purposes of tinder- 
standing them (called the semantics). For example , gome 
particular words (say, digits) are used to perform certain 
functions [28 ]. These words spoken hy a person could he 
stored suitably in a table or array in a computer memory 

and later on, when spoken again, they could be compared 

the 

with the elements of /table and if matched, certain chain 
of actions can be initiated. The linguistic aspect^ like, 
the stress on vowel, tone and intona'^ion are not utilized 
herein. However, it forms part of type (b). Though much 
work has been done in this aspect of speech recognition 
[1,4,9,21] yet because of the 'uncertain nature* (or 
dependency on number of factors like emotional state, 
sex, age etc.) of the speech, much more is req.uired to 
enhance the man-machine ccmimunication. Of iimaediate 
conceiTL to us here is the intonation. We would venture to 
extract the information about the intonation of speech 
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signal which can then ibe classified to indicate various 
types of statements. 

1.3 Definitions Intonati-cn 

Intonation nay be defined as the gross variation of 
the voice fundamental frequency over a sentence or phrase 
where a phrase or sentence is made up of one or more word 
strings in a given language. 

Both the phrase and sentence intonation carry important 
information from speech recognition point of view. What we 
shall be considering in the present work is the sentence 
intonation. As will be seen shortly the various intonation 
patterns would be obtained for different types of statements 
and hence their classification. 

1.4 Speech Production Mechanism 

The generation of intonation of an utterance is 
organizsed in parts, in terms of certain synchronized patterns 
of muscular activities of the larynx and the respiratory 
system. Intonation is therefore perceived in terms of complete 
contour of fundamental frequency and amplitude. The sounds 
are produced by a few basic mechanisms obtained in various 
ways. A regular periodic vibration is generated through the 
action of the vocal cords, which are .two bands of elastic 
tissues in the larynx. They may be opened to permit free 
breathing or brought together to produce various types of 
voice qualities. At a given moment, they are closed ti^tly 
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with narrow opening between them and if the air behind is 
moderately compressed by the contraction of thoraic cavity, 
the air pressure forces the vocal cords apart. A small 
amount of air passes through them and they again are closed 
because of their elastic nature. And thus the alternate 
closing and opening of the vocal cords give rise to voice 
fundamental frequency(f ^) . 

Based on the above facts the concept of source-filter 
model of speech production [7] was developed, as shown in 
figure 1. 

Accoustic Accoustic 

source filters 



e(t) f(t) 


Figure 1.1s A simplified model of speech production 
system. 

The vocal tract is excited by vibrations of the vocal 
cords modulating the air stream passing through the larynx. 
The variation in volume-velocity of the air stream is the 
quasi-per iodic pulse like waveform that is ultimately res- 
ponsible for the f\mdamental frequency of the voice. The 
vocal tract consists of oral, nasal and pharyngeal ■ cavities. 
It may be considered as a linear filter network with a 
fundamental frequency waveform as its input and the radiated 
speech signal as its output [8], 
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Prom the communication point of view [20] physio- 
logically in general, for males the f-undamental frequency 
varies from 70 Hz to 160 Hz; for females 140 Hz to 320 Hz, 
and upto 500 Hz for children. Health factors affecting the 
f^ here would he assumed to he inconsequential. The effect 
of emotions on the fundamental frequency f^ have been studied 
extensively by others [SO]. 

Sentence intonation is nothing hut the contour i of 
the fundamental frequency f^. Therefore measurement of f^ 
over the period of utterance would result in the formation 
of intonatuon patterns. 


1»5 Types of Intonation Patterns 

The various types of sentences may he distinguished 
from the study of the variation of f^ during the utterance. 
Terminal glide in the f^ gives an important clue to the 
type of the sentence. Questions are often distinguished 
from statements hy a teiminal rise in f^ as against a 
terminal fall in a statement. However, a question may 
also he distinguished hy a comparatively high f^ through- 
out the utterance P-6], Spectrographic analysis of speech 
also have shown that, questions tend to he spoken on a 
higher f^ than statements usually ending in a moderate 


rise P-6], These facts suggest that not only the direction 
and rahge of terminal glide (rise and fall) hut the shape 
and level of the entire contour affect th-O listner's 
judgement as to the type of sentences perceived. Broadly 
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speaking the sentences could he divided into the following 
mjor categories; 


1. horEial statement 

2. ..Questions 

3. Anger 

4. Delighted surprise (exclamation). 


As shown elsewhere [l2], the average, Patterns for the 
intonation of above 4 types could be illustrated as in 
Figure 1.2,. The utterance chosen was ’five thousand six 
hundr.ed ten' spoken in different ways. The variation of 
f^ at the terminal glide is marked in the case of question 


P 

TJ 


N 

W 

''55\ 

w 95 .1 
g 210-1 

S 1204 — 

^ 2501 

J-H 

byn 


120d 




SIA-TEMEHT 


QUESTION 






ING-ER 


\ 

SURPRISE 



200 400 600 800 1000 1200 

Time (m.sec.) 


1400 


Figure 1.2* Average fundamental frequency patterns of -Hae 

intonation contours# ,, , 



7 


Based on these patterns we can analyse the fre- 
quency component of an utterance with respect to time 
and depending upon the type we may classify the type of 
sentences. 

1 • 6 A Recognition System 

It will he appropriate here to see one example 
of Man-computer communication system through voice and 
orient our problem of intonation in this environment. 

Let us consider a block diagram of a general 
purpose man— computer communication system through voice 

[26] as in Figure l.J. 

It is clear from the block diagram that the input 
speech after digitization is stored in the computer 
memory. This is then analysed through suitable spectruim 
analyzer. The output of the analyzed speech is recognized 
as a string of sounds (phonemes) and then segmented to form 
sentences and suitable classification is done thereby in 
phase ¥; the message reconit ion interpreation of the speech 
is done and appropriate responses are generated suitably 
and outputed as desired bythe interac tor. Careful study of 
literature on speech analysis [ 8 ] shows that even with 
different types of spectriim analyzers often it is diffi- 
cult to resolve the ambiguities of recognition of a diffe- 
rent classes of speech sounds, e.g. a vowel /O/ and a glide 
/w/ might be confused on the basis of their gross spectral 
information - formants, duration etc. To resolve these 
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difficulties, one may think of two of the following 
approaches : 

(1) ask the speaker to repeatedly speak the message 
very slowly and clearly. However, this may not 
always he possible, 

(2) take the cues from other acoustic-linguistic aspects- 
say, intonation contours. 

The second seems to be viable because even though the 
recognition of a particular speech sound may be difficult in 
a given sample, since an intonation type is always present 
in a speech, this can be utilized suitably. To illustrate 
this point, let us take an example: 

Suppose a normal man with normal organs of speech 
(vocal set up) speaks a sentence in English: say (we have 
transcribed in normal orthography here). 

When did you come ? 

If on analysis, we note that the speech segment /w/ 
could not be clearly, and unambiguously (say undiff erentiable 
from /O/ or /v/) identified, what one can do is as follows: 

1. Look at the intonation contour. 

2. Does it raise at the end of the utterance ? if yes, and if 
the second segment is identified as /h/, and a vowel 
after it, then the utterance begins with wh - (who, 

when etc. ) . 

3‘ Otherwise call other recognition algorithms. 

Without elaborating further, we note that in English, 
normally a querry sentence ends with raising terminal 
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intonation, irrespective of whether a querry word such as 
what, when, which, who, occurs at the beginning or not. 

Thus one can construct a complex algorithm suitably incor- 
porating these facts. The important point to be noted here 
is, that the intonation analysis of a speech input could be 
really used to call appropriate types of spectrum analysis 
techniques (hardware or software wise) and thus enhance the 
recognition procedure enormously. Secondly, while the normal 
speech analysis techniques restrict their attention to the 
recognition of vowels, syllables /words, this approach woxild 
open procedures for continuous recognition of speech. To 
our knowledge, we are first to introduce this approach in 
speech recognition system. 

In Figure 1.4, we have provided the modified version of 
Figure 1.3 system incorporating the intonation analysis 
procedure before the spectrum analyzer is used. 

In this thesis an attempt has been made to establish 
the feasibility of the analysis, recognition and classifica- 
tion of intonation patterns as a first step in the direction 
of total continuous speech recognition. Using a bank of 
simulated digital filters for this purpose, a preliminary 
investigation has been done with the synthesised signal 
samples. In Chapters 3 and 4 digital filter simulation done 
on IBM 7044 is described and the problems encountered and 
the suggestions for further work are given in Chapter 
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REVIEW OP LITERATURE 

Speech research studies have gained momentum 
since I960. Comprehensive survey of literature on the 
automatic speech recognition system exists [l?]. This 
can he further supplemented suitably to include materials 
published upto 1973 (Appendix A). Most of the studies can 
be broadly classified as follows; 

(a) Isolated spoken vowel recognition 

(b) Syllable recognition 

(c) Humeral (digits) recognition 

(d) Word recognition 

(e) Speaker recognition 

(f) Prosodic feature recognition (stress and tone). 

Apart from the areas of speech recognition 
different techniques, hardware and software have emerged 
in hthe recent past. A few broad systems can be siommarized 
as follows; 

-(a) Hardware recognition system 

(b) Hardware exam computer combination 

(c) Adaptive (learning) recognizers 

(d) Usual accoustic analysis 

(e) Simxalat,ed recognizers. 
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In tlie following brief outline of the various 
recognition systems we shall mention only the most 
significe-nt ones and rest of them may be found in the 
litera,ture [ 17]. This however does not mean that we 
underestimate other approaches or overestimate the 
contribution of the systems to be mentioned. 

Before we deal with the various speech recognition 
systems, it will be in order to mention briefly what we 
mean by speech analysis which forms the preliminaries 
required to understand the ensuing survey. 

The usual way of visual portrayal of speech sounds 
is the use of sonagraph or spectograms. Spectograph (or 
Analyzers) produces a permanent spectogram of a complex ! 

signal. Developed at the Bell Telephone labs [15] in mid- 
forties, it is used in extracting many speech features such 
as formant structure, voicing, friction, nasality and pitch. 

2.1 Hardware Systems 

An example of the use of hardware system in the ^ 

recognition of conversational speech is given by Suzuki{5 0] -uBiJi 
a 26- channel spectrum analyzer covering frequency range from | 
200 Hz to 5»900 Hz. The channel outputs were separately 
grouped and connected to individual vowel decision circuits. 

The speech was segmented into voiced and unvoiced segments 
and envelope intensity measurements were also used. The 
final classification was made by observing the phoneme 
most frequently recognized , 
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A device for recognizing Japanese digits spoken 
in isolation [2l] was based on number of voiced interval 
in each utterance, the presence of an unvoiced consonant 
and formant frequencies and frequency slopes at three time 
points in each voiced interval. 

Another promising approach to real time speech 
recognition using analog threshold logic (ATI) gates has 
been developed and the RCA laboratories [22 ] . These gates 
give zero output until sum of excitatory and inhibitory 
inputs is less than a specified threshold. Once this 
condition is satisfied, the output is linearly proportional 
to the sum. The basic ATI gates are used to abstract 
relevant features from the output of the spectrum analyzer. 
Some of the various features extracted are formants, zeros 
positive and negative slopes etc. Such a feature detector, 
instead of a categorical decision as to whether a feature 
is present or absent, gives a q\ia.ntitative measure of the 
feature. Outputs of these feature detectors are combined 
in appropriate manner to perform recognition. 

A system for recognizing Italian numerals [11 ] uses 
a 17 channel spectrum analyzer followed by threshold 
detectors which reduce the range of outputs from each 
channel to binary statements. These signals are then 
sampled and connected to separate circuits for recognition 
of each digits , 
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2 . 2 Hardware- Computer Gombination 

Denes f 6 ] gives a system for recognizing the ten 
digits from complete word patterns. Outputs of a 17 channel 
spectrum analyzer were sampled and recorded on a magnetic 
tape. This is then inputted to a computer to form a time- 
frequency patterns for nuimber of utterances of each word and 
stored as reference. Cross-correlation process was used to 
compare the unknown utterance 'with the stored patterns. Best 
match pa.ttern was chosen as that of the unknown signal. 

Another system due to Gold [13 ] uses a 16 channel 
spectrum analyzer, a pitch extractor and a voicing detector. 
The computer program segmented the words into approximate 
phonemic units on the basis of voicing magnitude and spectral 
information. Further spectral durational and intensity 
measurements were then made on a group of five segments 
centered on one segment which was defined as stressed. A 
scoring method, based on the similarity with previous 
measurements on known words was used to classify the word. 

2.5 Adaptive Recognizers 

The system is trained first for an utterance by a 
single talker and then used for recognition purposes. A 
recognition system [ 5 ] using adaptive threshold elements 
incorporates a 15 channel spectrum analyzer, a 13 level 
amplitude quantization, and the output is sampled every ten 
m.sec. This is then fed to ten adaptive threshold elements 
simulated by computer. 
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Another system [ 9 ] divides the frequency range of 
500 Hz to 3000 Hz into ten channels, each of which is 
envelope detected. Two more signals are formed from the 
energy above 3 KHz and the overall spectrum. These 12 
signals are sampled for 4 * sec. 8 ,fter the onset of a word 
and the resulting 192 samples were punched on to a paper 
tape to be fed to computer subsequently. Word recognition 
is by means of linear decision hyperplanes which are adjusted 
to give optimum performance during a training session. 

Recently Clapper [ 2 ] has given a scheme for adaptive 
memory technique to adjust to each speaker. Separately 
spoken individual words can be automatically recognized using 
a two dimension pattern of spectral density vs, a nonlinear 
time base. The pattern for a given word differs from person 
to person and must be adaptively learned by the machine for 
each speaker. Simple circuitry is given that learns a word 
with a single utterance and recognizes it therea,fter. 

2.4 Usual Accoustic Analysis 

A highly successful recognition system based on 
direct analysis of speech wave has been developed by 
Reddy [ 27 ]. Direct processing of speech wave for segmen- 
tation and feature extraction is the significant aspect 
of this system. Speech wave after normalization is 
subjected to a segmentation procedure which identifies 
sustained and transitional segments. The main parameters 
used for this segmentation are intensity and zero crossing 
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rate of consequetive 10 msec, sections of the signal. 

Minimal segments (10 msec.) are grouped together if they 
are to be accoustically similar, similarity being defined 
on the basis that intensities of adjacent elements do.-s 
not vary by more than a specified tolerance interval. 

A vowel recognition system [30 ] performs the 
spectral analysis by means of an ’analogue ear*. The 
output of analogue ear was sampled at a high rate in 
4 msec, bursts synchronized with the onset of each glottal 
cycle.- Analysis was performed by a computer' by means of 
correlation operations. 

for 

Apart from these studies excepting/one study on 
intonation analysis [10 ] other studj^- of intonation pattern 
exists only from the psychological point of view. For 
example from the perceptual point of view the relation 
between the movement of fundamental frequency f^ with 
respect to time in giving rise to various intonation 
pattern has been brought out in [12]. This does not involve 
any recognition by a machine but throws much light on the 
perception of the intonation patterns. 

Thus, we realize the necessity or desirability of 
using intonational analysis through a computer system which 
does not exist elsewhere. ¥e describe our approach in this 
regard in the ensuing section ./ 



CHAPTER 5 


STATEMENT OF PROBLEM AKD PROPOSER SOTUTIOH 

It has been observed from the literature survey on' 
the speech recognition techniques that verj little signi- 
ficant work has been done in the field of intonation 
analysis. Total continuous speech carries important in- 
formation in the form of intonation contour. Hence the . 
need arises of the analysis, recognition and classification 
of intonation patterns which would be a step towards the 
enhancement of the work on speech recognition. 

The aim of this thesis is to obtain a procedure 
to be implemented on digital computer to analyse and 
recognize the intonation patterns of the specified 
sentences. 

The variation of the intonation contour is essen- 
tially a plot of f\indamental frequency f^ of the speech 
signal with respect to time, f^ varies during the utterance 
of the sentence, hence its , evaluation versus time by digital 
computation constitutes the crux of the problem. It has been 
mentioned in Chapter 1 that f^ lies in the range of 70 to 
500 Hz. for male, female and children speeches. However, 
for our purposes, we wo\ild limit this range from 70 to 300 Hz 
which includes the signals cf Bale and f emale speakers 



18 


In order to process tlie speech signal in the 
digital computer, it is essential to have the analog 
speech signal converted to digital signal. An A/D 
converter would give the corresponding digital signal 
for speech. The extraction of the information atout f^ 
at various time intervals can be done by having a bank 
of band pass filters, covering the entire r&inge of 70 to 
300 Hz. It would be desirable to have the bandwidth of 
filters as narrow as possible, hence giving rise to a 
large number of filters in the bank. Compromise between 
the bandwidth and the number of filters wo\ild have to be 
done depending upon the filters used. ¥e would be using 
the digital filters whose stability is a factor of the 
word length of the computer, the amoimt of round off error, 
and truncation error and the location of the poles of ■ 
transfer function with respect to the unit circle in 
z domain [25]. Therefore the n\amber of filters and the 
bandwidth chosen would be decided upon by: 

(a) Stability of the filters. 

^b) Time requirement for computation 
(c) Memory requirements of the program. 

The system which would be used for intonation 
analysis-recognition is shown in Figure 3»1*. The 
following assumptions are made in the analysis* 

(a) The .continuous signal to be analysed has the 
explicit boundaries {of the sentence) defined ' already. 
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(To) The duration of the sentence spoken is approximately 
of 2.4 seconds. 
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Figure 3*1: Proposed analysis recognition systen 
for intonation patterns. 

The system comprises of three main blocks! 

(a) Converter 

(b) Analyzer 

(c) Recognizer. 

The converter constitutes a transducer and an 
A/D converter. The transducer produces appropriate 
electrical signal corresponding to input speech, to be 
fed to the A/D converter. The A/D converter gives the 
digital output for the analog (speech) signal input. The 
speech signal would be sampled at chosen regular intervals 
and a binary coded signal would be available corresponding 
to these sampled values. The sampling should be done 
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at a rate twice the maximum frequency present in the 
speech. The output of the, converter is stored on a 
ma,gnetic tape to he processed on the computer later on. 

The blocks analyser and recognizer are simulated 
within the computer. The digitized speech is first passed 
through the low pass filter. ^A* to give signal sequence' 
containing frequencies below 500 Hz. This sequence is 
then passed through the band pass filters to BPP^ 

successively. The output of the filter bank is stored 
in the C matrix (code ma.trix), which represents the time- 
frequency relation of the speech signal. This is the 
intonation pattern of the utterance or the speech signal. 

A new set of computational procedure is required by 
which the contents of C-matrix can be classified by the 
recognizer. The standard intonation patterns for various 
types of sentences are stored in the computer memory. By 
selecting suitable algorithms the C ma.trix contents should 
be compared with the stored pattern. Decision on the type 
of the pattern should be made based on the above comparison. 
This aspect of the system requires a further investigation. 

However, in this thesis we are restricting ourselves 
only to the analyzer part and rest of the system is assumed. 
In the next chapter, we give the details of a simulated 
analyzer. ■ 
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SIMULATED ANALYZER , 

Processing of the digitized speech can best be 
analysed fT or various frequencies by passing the signal 
through a bank of digital filters simulated in the 
computer. What follows in this chapter is a discussion 
on realization forms of digital filter transfer function, 
the design criterion, simulation of filter bank and the 
choice of these filters for our purposes. 

The term digital filter refers to a system which 
executes an algorithm by which a sampled signal or 
sequence of numbers x(nT), acting as an input is trans- 
formed into a second sequence of numbers y(nT) termed 
the output signal. 

One such algorithm used very often is 
N M 

y(nT) = a.x(nT-iT) - b^y(nT-iT) (l) 

1^0 ^ i=i 

Knowing the present and (N-l) past values of input 
sequence x(nT) and the M past values of output sequence 
y(nf) the present value of the output can be computed. The 
difference equation (l) in time domain gives the transfer 
function of the digital filter in z domain as 
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^ .1 

2: ^1 

H(z) = (2) 

M 

1+ 'b-z’’^ 

where z**^ represents the unit delay operator. 

The transfer function H(z) given hy (2) leads to 
two types of filter realizations; (a) where at least one 
Uj^and one is nonzero - gives recursive digital filter, 
and (b) where all b^ are zero, which gives nonrecursive or 
transversal digital filter. 

Recursive digital filter output depends not only 
on the input sequence but also on the previous values of 
the output. Ifftiereas the transversal filter output depends 
entirely on the input sequence. Both these types have 
distinct characteristics. The transversal filter possesses 
excellent phase characteristics, but it requires very large 
number of terms (hence more number of multiplication and add 
operation to obtain a sharper cut off. However the recursive 
filter requires few terms to give a much sharper cut off. 
Thus to obtain same cut off recursive filter requires less 
computational efforts. 

The required transfer function H(z) for digital 
filter can be obtained in one of the i.two ways; (a) obtain 
the transfer function in s-domain for the corresponding 
analog filter and apply suitable transf ormation to obtain 
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z~domain transfer function, Ifiasodj ("b) obtain the 
transfer function of the digital filter directly in 
z-domain [25 ] • Both methods have their advantages. In the 
first case, the use could be made of the extensive work 
done in analog filter synthesis. Among the various choice 
of filters, a few are Butterworth, Chebychev and elliptic 
filters. All these filter characteristics are approxima- 
tions to a desired rectangular pass band. The Butterworth 
kind achieves this through monotonic amplitude versus 
frequency characteristics. By allowing a ripple in pass- 
band the Chechychev kind using the same number of poles 
and zeros can achieve sharper:-. cut off* Elliptic filters yield 
even sharper cutoff than Chebychev for same network complexity. 


Synthesis techn:^ r/ies for analog filters yield transfer 


fionction in s-domainj 


K 




H(s) 


k=o 


El d,s- 

llao ^ 


K <; L 


( 3 ) 


The transfer functionoof digital filter can be 
obtained from this by applying one of the following 
transformations: 

(a) Standard z-transformation 

(b) Bilinear z-transformation. 

In first case transfer function H(z) is evaluated 

log z 

using the transformation s-> — , where T is sampling 
frequency. 
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Tijus 


H(z) = H(s) 


iS= 


loggZ 


(4) 


This is conveniently done by breaking H(s) into 
partial fraction and tben making the substitution 


1 

s+a 


T 


1-z e 


Standard z transformation is good in preserving 
the impulse response but its frequency response suffers 
from the aliasing error. This error is due to the fact 
that H(s) is not band limited. However, if, it is band 
limited, i.e. H(s) = 0 for,: wj > Wg/2j w^ being angular 
sampling frequency, then aliasing error would be absent. 

This aliasing error, can however, be eliminated by 
using the bilinear z transformation or z-form. This 
requires the substitution 


1-z 


-1 


1+z 

This transformation, eliminating the aliasing 
error, carries a nonlinear warping of the frequency scale 
due to the fact that the whole of the frequency scale axis 
of s plane is compressed into a portion | | ^ 

can be compensated by prewarping the critical frequencies 
of the analog filter in such a way that application of 
Z“f orms: will, these frequencies back to their original 

values. 
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Thus 

D. 

= tan ( 5 ) 

gives the analog frequency for corresponding digital 
frequencies w-r, [ 25 ]. 

jL 

Another transformation matched z transformation 
could also be used for obtaining H(z) [l 4 ]« 

It is tLifeaal practice to design a low pass filter 
first and then using suitable trans format ions, high pass, 
band pass or band stop filters could be obtained. These 
transformations could be applied either in analog filter 
[31] or in digital filter [ 3 ] transfer function. 

It has been shown [3] that if H(z) is the transfer 
function of low pass filter, then H(l/f(z)), where 

^ \ 2 ( z ) ■ ' 

f(a) = ^ — sli-t 

. ( 1 -a z ) 

acquires a band pass characteristics. 

Here a is so chosen as to give appropriate cut off 
of the band pass filter, a is given by 

.cos k:(wp + w, ) 

a - cos { 2 k w ) = ( 6 ) ( 

^ cos k(w2 - w^j 

where w^ = centre frequency of band pass filter 
= upper cut off frequency coy. BP filter 
w^^ = lower cut off frequency hflEP filter 
k - fr ~ angular sampling frequency. 
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The cut off frequency for the corresponding low pass 
filter is given by 

W, = Wg - (7) 

Realization of Filters 

A recursive filter represented by (1) can he 
realized (or simulated on computer) using either direct, 
cascade or parallel form as shown in figure 4*1. 

The d irect form is the direct realization of the 
expression of transfer function hwhereas cascaded form 
is obtained by expressing transfer function as 

H(z) = H^( 2 ; X 112(2 X ... 

where each of the subfilters includes a subset of the 
poles and zeros of H(z"'^). Parallel form is oUttained by 
expressing H(z*’^) in its partial fractions. 

Direct form requires a larger amount of memory and 
greater accuracy in filter parameter determination. This 
is seldom preferred to the cascade or parallel realizations. 
However, cascade form is preferred over parallel form 
because the errors due to finite word length in computer 
is pronounced in parallel form. 

System Design 

The speech output of the A/D converter has been put 
on the magnetic tap Of for processing later on the computer. 



Direct^ form H(z ) =: 

l+F(z”^) 



Parallel form H(z“^) = H^(z~^3 + H 2 (z”^) + 


Figure 4.1 2 Various realizations of digital filter 
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The sampling rate of the signal is 20 KHz using a word 
length of 12 bits. Offset binary was used in coding. 

As has been pointed out earlier, the intonation 
contour is the variation of fundamental frequency ,;f^, 
of speech with time. It has also been shown that this 
fundamental frequency f^ of speech lies in the range of 
l0 to 500 Hz. 70 Hz to 160 Hz for male speech, upto 
300 Hz for female speech and upto 500 Hz for children. 
Therefore, in order to evaluate the intonation contour, 
it becomes necessary fhsttltiheesjgnal be passed through a 
low pass filter with cut off at 500 Hz. before processing. 
The resulting signal could then be analysed for intonation 
contours. Specifications for this filter would be 
Cut off frequency 500 Hz 

Pass band tolerance 25 percent 

Stop band tolerance 15 percent 

Transition ratio 0.8 

Sampling frequency 20 KHz. 

The type of filter chosen for this is elliptic 
filter. Since it would give the required filtering with 
minimum number of cascaded sections. A computer subroutine 
PROIP due to Sinha [ 30] was used to carry out this design. 
This subroutine returns the coefficients of the digital 
filter which is simulated as subprogram DIG-PIL to give the 
required output. 
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Filter Bank 

The fundamental frequency f^ of 70 to 300 Hz has 
been divided into contiguous channels of equal bandwidth 
of 25 Kz each giving the 10 channels of following bands. 


70 

- 95- : 

Hz 

195 - 

220 

Hz 

95 

- 120 

Hz 

220 - 

245 

Hz 

120 

- 145 

Hz 

245 - 

270 

Hz 

145 

- 170 

Hz 

270 - 

295 

Hz 

170 

- 195 

Hz 

295 - 

320 

Hz 


Here again the type of low pass filter used was 
elliptic filter and the subroutine PROIP was used to obtain 
low pass filter coefficients. A prototype low pass filter 
was designed first with following parameters: 


Cut off frequency 
Pass band tolerance 
Stop band tolerance 
Transition ratio 
Sampling frequency 


25 Hz 

20 percent 
10 percent 
0.9 
1 ZHz 


Sampling frequency in this has been chosen as 1 KHz as the 
maximum frequency present in the signal is 300 Hz. Achieving 
this, however, is no problem since 1 in 20 samples of the 
input signal can be taken to suit this condition. Transfer 
function of the prototype low pass filter is obtained as 




n 




where n is the number of cascaded section. 


( 8 ) 
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For further discussions let us take n = 1, which 
however can he extended to cases n 1 without difficulty. 

Converting this low pass filter to hand pass filter 

hy the transformation z‘~'^ ^ in (8) 

(1 - az"^) 


.-1. (1 - az (1- az^) 

H{z ■^) = e-, IT— ZT I1- - I 1 ^ (9) 

(1- a z (1 - az 


Upon simplification, one gets 

-1 -2 -3 -4 

, Ctt+CtoZ +CttZ +Ct.Z -^+CTr-z ■ 

tt ('_-*1 n 11 12 13 14 15 

tnz j = 


( 10 ) 


where the coefficients in (10) are related to that of (8) 
hy following relations 


11 ^ 

^11 

dll 

12 

a(a^2~2a^l) 

^^12 

13 

a '• (a^i~ai2+®-i3 )“ai2 

di3 

14 ” 

a(ai2~2aij) 

di4 

15 ^ 

^13 

'^15 


a(bi2-2bj^^) 


'11"“12'^"13' “12 


(X(bi2-2bi3) 


( 11 ) 


Thus one comes to conclusion after observing (11) 
that the coefficients of various hand pass filters can he 
obtained from a low pass prototype filter coefficients, 
the only parameter changing is a, which is decided upon 
by the cut off frequencies of the corresponding filters. 
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IJransfer function in (10) can be represented 
pictorially as shown below. 



This has been simulated on digital computer using 
the subprogram DUFIl. The tinit delay z ^ represents 
storage for one data sample. 

Testing of Pilters 

Before using the filters for actual processing, 
they are tested for (a) Impulse response (b) frequency 
response. 

ImpiAlse response indicates the stability of the 
filter. This could be obtained by getting the output of the 
filter for an input sequence of 1>6,0 ... 0. for a stable 
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filter the output decays in transient state and becomes 
zero in steady state. 


Frequency response of the filter is obtained by 
1 • j w i ’ 

replacing z by e in the expression for transfer 
function, i.e. from (10) 




c +c 
11 °12^ 


^11**^^! 2 ^ 


;+Ci^ e " e e 


■+d^^ e 


+‘=14 ® 


■+di5 e 


( 12 ) 


dividing numerator and denominator by and simpli- 

fying gives 

C (c^2.+‘=15 

C (dii+di^ ) cos 2 wT+ (di 2 +di^ ) coswT+di 3 ]l 

i-3[ (ci2-Ciq)sin2wT+(c i2-Ci^)sinwT ] 
+3C (di2-di^)sin2wT+(di2-di^)sinwa? ] 


H(w) for various values of w can be plotted to find whether 1 
the response is the desired one. 

Program Imnlementation f 

[ 

The aperiodic nature of the speech wave necessitates ■ 
the creation of a data window. This subset of data (i.e. a 
portion of speech signal) is passed through the bank filters 
respectively number of times till the filter output stabilises j 
It has been found by experim^tation that 5 such repetitions | 
give the desired results. The data window chosen was of 
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50 samples. Once the corresponding output is obtained 
the r.m.s. value of these 50 samples is taken and compared 
with the r.m.s. value of the corresponding input samples. 

To ascertain whether the output of a particular filter 
is present or not, use is made of the fact that filter would 

attenuate the signals outside its pass band* Whereas the 
signals in the pass band would give larger outputs* Thus if 
the ratio of the r.m.s. values of the input and output of the 
filter is above certain threshold then the presence of the 
particular frequency is indicated in that particular subset 
of speech signal. 

Thus the decision taken at the end of the processing 
of 50 samples, based on the value ratio will be indicated as 
1 or 0. If the signal is above the threshold, indicating the 
presence of frequency in that period, would be indicated by 
1, and the absence by 0. 

The next 50 samples are taken subsequently and 
processed similarly, till all the samples are processed. 

The signal is then passed through the subsequent filters* 

The decisions 1 dr 0 at various time instants are 
stored in o-matrix. 

The computer program organization and the flow 
charts for the above program are given in Appendix B. 



CHAPTER 5 


RE SUETS OE EIPERIMEHTS 


The low pass filter A (Figure 3 .I) for speech 
signal and the hank of filters were designed and simulated 
on the computer. The coefficients (i.e. a^’s and hj^’s) 
of the low pass filter A and that of prototype low pass 
filter are shown in Table 5.1* The coefficients of the 
various hand pass filters were obtained from the low pass 
prototype filter using equation (11 ) of Chapter 4« These 
coefficients (c^'s and d^'s) are tabulated in Table 5*2. 

The testing of the simulated, analyzer was carried 
out in the following manner. 

A signal comprising of five sinusoids of different . 
frequencies was fed to the program package. Five filters 
in whose pass bands these sinusoids are covered indicate 
their presence by 1 in the relevant row. Whereas the 
absence is indicated by a 0 in the row for which no _ 
sinusoid existed. The results are shown in R.l. Another 
similar signal, but the frequencies of sinusoids staggered 
from that of the first signal was also fed. This results for 
this are shown in R. 2. These results establish the efficacy 
of the filter bank. . • 

Because of the unavailability of the digitized 

signals were synthesised by mixing several 


P 
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sinusoids together so as to resemble the spectrum of 
the speech signal. Sinm.oids of different frequencies 
appeared at various tine intervals in such a signal. 

Two such signals were synthesised and were analyzed by 
the program successively. Results for these are shown 
in R.3 and R.4- displaying the c matrix. They indicate 
the presence of frequencies in the various time intervals 
and in corresponding bands as per the synthesised signal. 

It therefore can be inferred that similar type 
of results would be obtained when a digitized speech 
signal is given as input to this program package. 
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CHAPTER 6 


COHClUSION 


Th.e intonation contour is an impor'^^ant feature 
of speech. Its evaluation would prove to be a signifi- 
cant tool in the continuous speech recognition systems. 
In the present work a procedure has been developed to 
find out the intonation contour of an utterance and the 
feasibility has been established. Some of the specific 
contributions of the thesis are: 

(a) Simulation of bank of digital filters. 

(b) Use of frequency transformation techniques 

to obtain appropriate band pass filters 
from low pass digital filter transfer 
functions » for the purpose of speech 
spectrum analysis. 

(c) A procedure to evaluate the intonation 
contour of a continuous speech signal. 

The program package developed, gave satisfactory 
results with the synthesised speech signals. However, 
it could not be tested for actual speech signal. The 
magnetic tape containing the output of A/D converter 
(by courtsey of T.I.F.E^Bombay) could not be used 
probably because of the incompatibility of the systems. 
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Attempts were made to design band pass filters 
with a narrow band width (10 Hz). This, however, resulted 
in oscillations at the filter output. It has been found 
that at narrow bands, the poles of the transfer function 
become very close to the unit circle. Due to rounding off 
error in computation, it is possible that the poles fall 
within the unit circle, leading to instability. This 
difficulty was, however, overcome by choosing a BW of 
25 Hz. .... 

Scope for further work 

In order to obtain better intonation patterns, 
narrow band filters should be used. It is therefore 
desirable to obtain stable narrow band filters. This 
particiuLar aspect requires further investigation. 

Another aspect which requires further investi- 
gation is the selection of suitable algorithms for 
matching the stored patterns with that of the analysed 
ones. Once this has been done, the classification of 
the utterance on the basis of the intonation patterns 
could be achieved. 


Speech Recognition and Analysis Review 


o 

m 

H 


Q) 

U 

o 

0 ) 


0 

S 

O 


k 

O 


Q 

Ph 


m 

0 

o 

« 

0 

0 

<M, 

0 


0 

0 

iXi 


Jh 

-P 

o 

0 

A 

a:i 


!i 5 


0 

mD 

0 

P. 

Oi) 

d 

0 


0 


w 


<H 

O 


k 

o 

•rl 

-P 

‘H 

m 

o 

o 

0 


0 ♦ 

.H 

o rt 




d 

s 




■2 

.3 

Hi 


urv 

V0 

':Cr> 

H 


4K 

mmunication 
, 11.4 


O CM 

O 

O • 

-H 

H 

-P 

Xi H 

P4 

5 ::^ • 

0 

cd 0 

O 

0 

u 

0 03 

0 

0 

Ph 

•H 

nd 

0 H 


5h H 

05 


0 

0 ^ 

X^ 00 

•H 

>>CM 

0 

0 

0 

m 

rd • 


4> lA 

•H pH 

ri • 

^ • 

&“" 

p-[ ^ 

0 ti 

P^ 

d 0 0 

0 o5 

S) 0 -H 

*H ;cl 

•H ‘H -P 

0 O 

03 > cd 

>> 

P H 

H 

X{ 0 0 


^ CO 

d cd 

0*0 


H 

P 

Xi 

o!> 

0 


tr\ 

ft, 

OJ 


0 0 P 
p^ 0 EH 
ca pci^ 


0 p 


0) 

,CJ pi cv 

Q<a> 
05 fciO 



--P H 

•v ‘p 

• ri 

• k 

«! -H -d 


• 0 

*00 

S cd 

% 

.•k 


I. L 

r^ENlJi^AL 'L 



42 


0 


p d 

o 

43 

d o 

o 


0 P 

CO 

0 

o -p 

« 

o 

d 

d 


0 d 


0 

P^P 

o 

p 

P O' 

o 


O 3^ o 

0 


OP 0 

rt 

•H 

HP d 



d -xi 





a © 


d 


0 

;b4 © 

0 

0 


d 

© d 

0 

0 


o 

© 

fd 3: 


•H 

Pi© 

© P 


43 

© p 


0 


0 

•H 


,Q 

m 

H 

0 ‘*0H 


U) 

d 

H H 

d 

0 

p 

43 

t0P 

cj 

0 

•H 

0 

d 

CD 

P 

tcD 

0 

H 0 H 

VU 

•H 


1 OH 

0 

HP 



0 

P 

p p 

d d d 
0H 0 

0 

0 

d 

d 0 

0 

0H P 

H 

H0 d 


0 

CTv r 


cn p 

cr^ > 0 


© d! 



© P © 

d -H 0 



(D l5DO 

>1 -rt |6S 
© <© 

© © 



0(0 fH 

©Ho 


0 

a 


d 

0 >1 


0 

H H P 



ip 

on 

no 

'« 

H 

d 

p 


p 

H P P 

d 


p p d 

cJ 

H 

S P 0 H 

CM 


, 0 

0 

p 0 

0 • 

0 -H 

d d 

© © +3 

d 0 

0 p -H 

a i> 

d.(Q p 

d H 

© © caD 

0 m 

PiH 0 


H 0 

d P 

00 © 

0 0 

cr> © p 

H d 



0 


0 1 


HI d d 


P 0 p 

0 

d ^ d 

0 

H d 0 

H 

H 0 d 

P 

H H 

d 

0 0 

H 

0 H d 

H 

d 0 p 0 

>> 

0 H H H 

0 

a c0 0 43 


d d d 

0 

0 -H d‘H 

0 

H 0 0 0 

H 


•d 

o 

a 

-p 

© 



CQ -P - r- ^ 

© p> H Td © 

O -H a 
a tio M tio 

■H © P> © 

cs Td ra d © 

0 Q) 

ra © an 

•H T3 H © © 

m ^ h 9s 

ra o cd o p o 

<jj 3: ^ tiD ra > 


0 

• 

0 

0 P 

P 

0 0 

d P 0 

0 

OH © 

CH 

© © 1 > 
0 0 -H 

0 

d H P 

> 

P10 © 

*ri 

P 4 > 0 

P 

© 

0 

H d 

d 

!»©^ 

H 

d ©^ 

43 

© T-t P< 

0 

© 0 0 

H 

•H t> P 

pt piw © 


fcjD © 

m d Td 

d H d H -H o 

© H -H cd o 

Pi © d ° ® 

LPi o H S o m -rf 

© -H Pi-H © 

•>H LP, © © © 


© rH 


d fp 

o d 


ra cd N d 

Td^OOOPdH-P 
d d— '-P d! d © © 

3 © CO 4S’ 
d* P“ © 

. d H'-'Sh 
O^.HP W-H 
© © -H 

d ©H dH d-P © © 
SHan-PH® ©d! 

■H> 5 ©©dOd-P 

OOP>©ao®-P*lOO 


O O 
O 

-P d 
h 05 ' © 

& 3 d! -P 


P> 

© d 
H © 

PiO 
LTV a-H 
03 ©Sp 

•'H ©-H © di 
ra cd d ft-H 
•td l> © fciO © d 
d d a-H _-P 
cd © -H © d © 
d -P -P 
• d 'H 
O^-H P-H 
© R 
d © © >a 

p a ®H ~ . 

•H d d rd o o 
oo p H o o Pi d 


H d 
cd o 
d ’H 
p d> 

o d 


•H P • 
Td 'd 
© © 
fct0 d o 
d © -H 
cd ^ P 


© 

P 

d 

© 


di 

Pi 

21 

A 


4;»' 

H f 

t«b 

*H 0) 

rt d ^ 
0 ^ ^ 
M O *p 

o o u 

A 0 CS 

CQ « W 


nd 

0 

t 

O’ 

S 

H' 


i 0 
!> 

•rl 

'o d-S 

0 0 -P 
0 0 

-T-t 5* 

ca a ^ rO o 



0l 

4^; 

m 

4^1 

P 

0 

oj 

ml 


m 

^ 

© diM 

0" 0' P H 
*H 'OS 

> H «P 0 
'0 0 -H 
ft 


qd 

i © d 

© ^ d-^ 

^ *p 0 • 

0 p p -P !> 
p d o 0P 
-rf P ^ 0 P 


o (d 0 
0 H P 
H 0 0 
PfflP 


I 

0 

Td 

1 


0 

> 

O 



I 

0 

Id 

1 


cdf 

0 

H 


CM 

tf\ 


'0r\ 

'H 


■'0£>' 

'' '!£>" 

, r^: 


VO 

'''irv 

<3r\ 


VO 

m 

ch 

■H' 


r~| 

VO 

Ch 

H 



0 >s 



-H 0 P 



H -^P Xi 



HO)© PIP • 



Cd 0 -H 0 OP 

-P 


Jh Pi pi -H tJ -H Pl 

P 

0 

0 3 fciDP 0 P 0 

0 

tS 3 

> a S=I C 5 CQ -H 0 

0 

•H 

0 Ph -H a cd pl Pl 

P 

P 

0 H M © t(D 0 

0 rH 

tiD 

Pi P 0 Pi 0 ft 

Pi 0 

0 

0 Pl P P 0 0 


0 

O ' 0 pS Pl Pl © 0 

ir\ 0 

0 

ft QiJii -H -H Pi CM 

<T\ i> 

P 


rtd 0 

Pl *• 

P 


P -P 

•H H ora 

0 


cd p 

W © Td Pl -P 


cd 

© 0 Pl sd sd H 


CQ $4 

H 0 cd^ 0 

Cd 


H 0 

0 5:4 ^ 0 a p> 


0 m 

? 0 H © © -P 


1 

Is P 

0 0 cd © H © 

Cd 


0 0 
t> 0 




0 H a 0 red 

0 



H cd 


f 



10 •> 

P 



© -H -P 1 >> 2 

0 H 



-POed *s:l^ “^a 

cd 



b©a ••>iOcd H© 

p p 

f 


ftH -H H -P -H H P> fl © 

'H -P 

p> 0 poa 


P> + 30 -H +=3 flo© 

0 

p 0 0 P 

rd 0 

!dH©Om-Hfci 0 ®ffi;sd a 

^ 0 

P+^ 43 

0 rM 

0 ®fHl-HplPl©rSlT 3 ft© 

0 Pi 0 

-P 0 

P P 

is ‘©©EiDlH+^M 5 d 

P 0 Xj 

P reJ cd pP 

p Cd 

HO ©>+30 *+30 

0 0 

0 a Cd 0 

0 D 

o{>P'idPio+=ftHd +4 

-P P -P 

0 0 P 

cd P 

id •HpJ-H©©o 5 ©ft 

CQ 0 cd 

p 0 0 r-f 0 

0 0 

sd 0 u g b 

^ a 

0 H P 'H 

a-p 

cd+^OTd© ©©0 Hp 

0 Is 

S H P Cd P> 

fH 

^ ?H©ft 0 S: 3 H©Q© 

^ p 

CD a «N cd 

l>s*H 

Op-HOOpOed f+o 

p> p p 

5h [> cd 00 a 


OO.HHCd-Ht>H'Ha 

0 0 

0 p 0 -p p p 

Cd 

ft© o©i+mfH© 

nd ‘H -p 

-p 0 0 

0 a 

dSs:lpll>l>cd-H®opcM 

0 43 43 

H Td n:d 0 

p 0 

o<do ~idPiO+^-H© 

0 cd cd 

•H 0 0 0 0 p 

H P 

JnOHOOOOflOOf-l 

Dap 

p -H -P p -H 

cd^ 

djoraft raH-Hf>ao 


-P +3 CQ 

a d 

(]> Q jZj 

o - o cd 

!-h ra Pi 

CD H CD O 
P 4 0 ft CQ 

is 

O oo o 

(Tv !>t>- o 


a 

0 

s -p 
o o 

r=i 

p^ 

-p 

a> :s 


-p 

o 

a 

CQ 

0 H" 
|5 "H 

o cd 


a ftH-p o 

O p H 0 Pi 
02 -H CdPl^ 


O 0 0 

•H pj 

•P tjo pi 

r 5 0 iM • 0 

p cd 0 cC 

P! © pl 0 0 

0 a 0 fn a 'H 

a -p 

© S!J 3 p 

© Pl Pl cd 


o 0 Pl 

Pl - 

■H i>s is 0 ' 
cd P P P 

•H -H +> -rl © +s 

£)fl! 3 ©©^H,Ocd» 
Pl O © Pl O 'H Pt^ci 
•HHOO .^0 © 

■ Pl -H © Cd o HP 


^ o 
© © 


o p p P p! (S cd 


cd 

O H 

u m 


d © isjp P) 
H - P -H 0 
bo © tiDpl 


© 4 :|©cd 0 pT-l-H® 
tbp Nf> PiO'cJ'd^ 


f-t S 

0 0 
tsj 43 

>>-P 


o 

'H 

a I 

0 -H 


H 'x:J 

oj Ph ^ O 

« cd rCt rO 

Cd H 9^0 

cd ^ P 
^ JL, 0 rd 0 
O 4^ ^ 0 0 

0 O O Jh p -H 
0 0 -P O -H "P 
pt( p 4 cd *P cd ‘P 
mm B m 9iH 


M 

H 

® cd 

a 

si 

cd P 
^ 0 
o -p 

iH 
IjCS -H 
to CH 


ri 

cd p 
0 

0 ;:d 
CSI Pi 

>:»a 

fH O 
Cd O 

p 

cd cd 


©H 



a -p Pl 


1 1 

VX) p P -H cd 


•H cd 

cvj Pi p PH -0 


pi 43 

P ftH 0 0 0 H 

54 

p p 

©OPPIJhHP© 

Cd 

p 0 

ft © b cd -H -H cd © 


•H a 

>1 ft 0 0 is H b 

P 

0 t !0 

E+ © •> fH 0 

0 

0 

Pi © Pi 0 Pi © 

IS 3 

0 0 

Pi H © -H ©P H 

•H -P 

© © ©ftHlcdp © 

« 

cd fH 

*0 g >s'H P ^ cd Pi 

m 

° ^ 

oSH-Po©ad 

0 

edft Pl 

oededooftPiro 

0 

ft 0 

oH 3 ©a©o^ 

0 

0 ©H 

!> ©■^Pimra'HO 


m pp 


pq 0 
H 
H 
‘ O 

P 

m o 

0 - H 
P' p.p 

0 O 


§ 


m 

H 


U 

0 

43 

P 

P 

a 

o 

o 


P 

cd| 


0 

H 

d 



0 


0 

Cd 

•H14 


^ rd 0 

P 

j 4 cd EH 

cd 

(d D 

H 


000 

«H cd cd 

•H *H 0 



;4 cd • 

P P 'H 

N 0 

0 0 D 

^ cd 0 

fXi ^ 

m^p:^ 

CJ> 

jH 

ir\ 

VO 

<T\ 

a>i 


rH 



> 

•H 


1 P 

1 pJ 
d P 

43 0) 
H -H -P 
d ^ o 
4 0 >3 


C\! 

VO 

cn 

H 



•p +3 

jzj ^ S4 

<0 f-t <D U 

o <r> o 0 

-H -H 
0 'Pi 0 «H 

P<-p.H ft -H 
•H -P -P -P • 

0 ?50 S 4 cr\-H 

• -H 0 • tiD 0 
CTi -H tJ 

Oin3 -rl <T\ 


W 

tl 0 .H 

O 

■H fciO 
m 

iH 

9 H 
•H -H 
fn 4 ^ 
P^ CQ 


CQ 
Q 
1 O 

• cd Sh 
W CD K3 Cl> 

a 

5h 0) <D 03 
^4^ -P m O 

jLj H 03 0 ) -p'— 

an ^ vjd 
03 CQ 
^ ao 4 ^---s 
O O ‘H 

•H <^C> 


(D 


o cd 

CH 43 

T 5 CQ 
0 0 


N 

'H 

Ih 

43 

05 O 

o ao'^ • 

cvJ urs^ X 
-p cd*H 

•v 0 •v ov fn 

m 0 m cQ 4 ^ 

0 *r~» Q CD a3 

Or^H O a 

S 4 Oj 03 

o 5 a cd-p 


(an cqq ^ m 

0 o 3 *HC 3 D 

a o Jh 0 "P H -p 

O 0 0 -H 0 P p.}P o3 O P ^ 

^< Hramj 4 0 0 ?: 5 aH;=^a 3 


0 

o 


0 U 1 
OP 

^ s 

0 o 3 

m ^ 


I a 

t 0 *H p O 

JZ{ CQP CQ 

H 0 g $:! 0 

IS3 P^ o p 0 

H -H P P • t 0 P 


O 

•H 

P P 
•H O 

tiD U 
O 0 


H ^ 0 O 
O O 0 


03 P 4 S 4 
0 O 
rcJ o 

P cd P 
0 

000 


Xh ♦ 

000 

P 41 ? :s 
o3 o 0 o 

03 p 0 

0 0 000 P S:0H 0 a *H 

0 <zj o 3 OP 

•CP0VD 0 

H H P 0 m 0 

O 
0 03 


0 
P 

O O 

<=! H cdP 0 
0 0 0 > o 0 


O P H o 


0 

U 

0 

M 


cd 

S 4 ^ g 

© ® p 
© 4^+5 
^ P -P 'H 
-P "H O 


® 

tsl 


P . -P ^ © © >3^4 © O 
cy'cr's 4 ©-pf 4 ^+j -p 
© © -H P> © P 
fH U o p'rJP>'P'H-P'p 
tH 'H ft-H ce © o rs © 

© © 40 ft +3 


W 

?4 • © 

ra © 

f4 'H -P fH 
O -H 

tiD -P 
ca -H 
!5 ri 

O P 

Jh ?:! 

0 0 


PP 0 p}HPo 30 PO><}P 


rfo 5 a 0 pj O}Z|pJ 0 P 
a a -H o p ^ p o o 

' }h P P o 3 O p 


00 O 0 G 3 rtPP 

fx4pi^>ppqp cd O 


JZj Jh CD p 0 

3 P P H riij 

o cd cd o p 

o a o > H 


0 
0 
0 

cd *H 
P* tiD O 
cd -H 0 
hD P CEl 


o 



P 'H r 9 

■H'w. 0 

• I 

t>>P 

cd 0 
^ 9 P 
0 Jh o 
Jh o p 
a 0 o Td 0 p 0 


0 o o o •- 

p CD 0 

0 4J> CD p 


P 

tiD 

P 

p 



tiO >> 

cd rH 


0 P O 
0 P ^ O 
p P 0 0 P 0 P 
0 P Ih ?4 ;d 

^ p a 0 H H 

>> p--^P P cd P 

H O ^ 

O a p P E 4 O rt 

Jh o P H o • P 

pS fd o cd^P 

0 P pjzj J 4 0 0 p 

JSip 3 rl 0 PPOP 

pp 0 p ^ p ;:{ 

OP 

0 P 4:! 0 

poo P'--' 

I 

o 

p "H p P ^ 0 


Cd si 

0 p 


cd P 


t> 0 a 

o • 

O 0 p 0 

0 0 P O 

S-i P p 
0P 

. . . ^ P S cd 

P O p 0 cd cd fn 

Si Sh ri 0 ^ ^ 

H p p [3 P P 


cd 


P 

0 


§1 

cd ?4 


c\i 

I I 

cd Cf o 
s:! P 
cd p P 

43 •» 
H JZl 


r 


£500 
o p 


p 
p 

•H 

c3 O • 4 ^ 

cd iq-H Jh-H -H 
<a} O 

O O 4^ CQ ?4 

.„ . . H Cd © -P EH ^ r! -H -H o m 
rt+> y^ 4i !3DH Cd<i;cC!©CDpi4J-P 
OOffiO-,C!f 4 ©aSO-HpS 

©c: 3 s:^Ocd^Oi:^®cdf^oft 

C3^ftt>,pj3-P'HO,i4rH fs-rH Nri 
rH © iH ^ © !>■ £4 © 4 ^ O © -H 


a -p 

© 

© B 


15 O H> 


Pi O 
Cd 00 
+> H ft • 
cd cd ft o 
tiD -H © 
cd -P S H 
^ ©'-.pq 


rQ 

Pi Cd 

-H H 

43 p 4 

^ cd <i< 
cd O 

S-P Pi 
0 '--' 


UD 

Ch 

p 


MD 

(T\ 

P 







P 




a 0 


P 




ra a • 




p 

P 

43 43 © 


P 


d 

d 

43 ra • > 


P 


0 

0 

d © !>jt3 -p 




0 

0 

,© P ra © M 


P 


u 0 

d • 

d p 


d 

0 

0 p 

0 

t^O © d P 43 


0 1 

d^ 

p^ 

ft 

d !> 0 > ,p 


0 0 

Saw 



p d 0 d 


fH P 

p 0 

CM 

CTi 

ra -H © d CQ 


0 1 

0 d 

CD 

<y\ 

■P d ft © 


ft 

Pd 



a ra © iC4 



P d 



-© 34343 


0 

0 , — ( 



iM P © 0 d 


0 0 

ft ft 



fei 43 d d 0 


p Ch 

' — ^ 


0 

ra 





P 

>3 1 © 

>:► 



4K 

0 

PH d ra 

5^ 



0 

0 ^ 

>3 d H 1>3 P d H 

d 


0 m u i 

0 

p 0 

H © d H OP 

d‘ 


0 0 P 

0 

0 > 

© d d rd H 

-H 

0 

0 JLj ?Li P 5 ::^ 


§ 


CQ 


„ c\J I o 
Ss CM -P {> ^ ^ 


OJ •> i:5 I 

+3 CQ CD CD 

+= <i> Ih a OQ tn 

pS iH CD 0 CD 

0} ch iq CD 

p a 'H o o ^ 

fc3* 'H ^ rH ‘H 

OMD T3 p< ro H 


-P 

t 3 03 

G I ^ 

o3 O -P 

CQ ta C 3 

CQ -p ra 

C o cd 

o 2 -H ^ iM 

•H a -P O CD 

+* s:! cd +3 

o tiD'H H 43 
cd -H ?H cd 03 

OQ 03 o P 4 

43 to |> -H 
CQ cd (Q 'H 

^ 43 >s o 

cd 0) 43 ^ 

M qd ft d 
CD cd -H o • 

d -tf H CD -H d 

d -fi CQ ^ CQ CD 

43 tio +> -H ch 
cd -H d d 'H 

CD CD o d CD -H 

ft > 'H T-i 43 tJ 


43<ri •H43cddH-c2 

CD©43cd-HOO03 
H430S;2-H.H-Pd 
ft d d ^ -f 5^ * 

a ® • d cd cd © a 

oaa'dda5ioa43d 
CD©0©0©-HcddO 
H 43 43 o ;d 43 -d ftCH 
43ftCQtQC|43rad-l^ • 
" ' ■■ o cd © 


o a >j© 

{s; -p ra 43 o 


a 


d 

cd 


t>5 d 
CD 

© N 

rd d 

43 © H 
d cd 
d -H d 
o cp cd 

•H © 

43 d 
d 

o ©^ 
ft o o 
cd 

ra fttp 
© ra o 
d 

cd 43 © 
H d43 

ft © d 
I a ft 

d © 43 
© d d 
ft d o 
>> ra 
trj cd ra 
ra xi 
u) 943 


© 

XI 

43 


I 

•H 
* 43 
o d 


© cp 

, . . a o 

d cd d 
© d 43 m 
C 3 cd d 
© © 
43 


43 cp 
•H O 

d 

EiDd 043'. 
d o d d d 
■H fc(D-H ft d d 43 O 
a c«3 043 ^ Cd cd 
a -H d 43 d >5 d 
d d d o -H «p cd 

CQ 43 © 01 X -H X 

ClD ® © 

>3^3 !3D.d CQMD 




■H 43 
i d 

d 


d 

© 


© m © 
43 d *4:3 
d ra ra ^ 
© -H 43 
ad -H 
©•H d 

HH O 
ft ft d © 

a a -H ^ 

M © O O 


CQ 

d d 
ft-H 
, .d 

d 43 o 
© d -H 
P o ^ 
© d 
tiD-H H 


8 ,“ 

H TzJ M 
Wt ^ O 

05 

m 
j:? 

o ^ 

. rl ^ 
^ *P O 

k 5 


( — 1 05 CQ ‘H 


ft 


Q> U 

ri *H 05 rCj 

ft -P > -P 


H 

d 

u 


0 O 

a dP 

ps a> p 
^ ft?^ 

^ CQ 

15 

000 

fp H ^ 


0 0 W 

^ P O 

d P o 

o 0 o 

ft ft O 


m ?:::! 

0 0 p 
0 p 

a ? 

d ^ 


a 

o 

H tS 
0 t^DP 


0 o d p 


1 


1 

P 0 


d 

0 P 


P 0 

d 

0 

d P 

p p 

p 

• 0 p p 

d CQ 


0 ^p d 

ft 

0 -p 

P p d ft 


p 

o 


>5 0 0 

Jh 0 0 

d ^ 

P 


a 


0 

^ 0C> ^ 
o -f o 

0 *r| 


0 p 


0 

P 0 0 


p 


;d P 

0 

0 

0 d 

P 

'H 


p 

P 


0 0 

0 

0 


0 

ftp 

d p rd 


d 


0 p 

a 

d 

0 0P 

0 

o3 


0 

CQ 

d 

•H 'H 

0 

d 


0 0 



43 0 H 

d p 

0 

d P 


0 

•H CD 0 

0 


d 

p 0 

• 

P 

d nd 

d P 

d 

0 d 

0 

d 

W) CQ 

d 

d 

0 

^ d 

p 

d 

0 © -H 

d 

cd p 

d 

d 

qD 

0 >3 CQ 

0 


P 

fd rd 

0 

P 

© © © 

0 

m 

cd 

EH 0 

a 

0 

d S ,o 

0 

d 

ft 


t*-N 


H 

0 , 
d a <i> 

a 3 N i 

d U >5C<^^ 
^ P P *0 


P 

0 

d 


• 0:5 0 
PP ^ 
000 
« .d p 


^ p >5 0 o 
o o p 0 0 



© d © d 

© © d 43 PI 


CM PL|<il g © 

tr~ p d rd ra EH 


H ra d 

• 

H ra © EH P 0 

d 

• 1 — 1 


0 

cd > P . 

P 

P 

PPM 

P 

P 

iS3 d 

P P 0 

d 

cd 'p P 

p d 0 

ft 

p 0 

P cd S 

in 

VD 


UD 

MD 

KO\D 

0 ^ 

cr\ 

Chcr\ 

P 

p 

p' H 



46 




p 


0 


0 


P 


0 


Pi 


in 



0 


p 


0 nd 

• 

H 0 P 


&iD4^ P 


Q) 


U 

O 


c| <ri Cf-( 

4^ :p H -H Q *0 

a ^ ftco (a OJ 

O 02 -H 03 -P • 

9 H ra O rs cj 

•H O CQ Pi O S 

Pi-H 03 03 

03 -P a -P 03 +3 > 5 +j Pi 02 P> ^ fH 

ftrifnCQa-PlHPl-POpS-POO 
Woa3-H03c6O 03a3-p^ c3c3P> 
m o ptiTl 02 Pils a P( 02 eH p,® 02 


C 3 

. ri -p 

® a a O CQ p 
Pi -H d -H P O • 
Pi-P Pi -P '03 Pi Is -p 
o ® ® ® 


Pi 

•H 


I OQH 

I a 02 r— I ® rH 

Oj S rH ® ^ ® I — I 

Pi a 03 Pi ® P H 

P PH Pi 03 

>3 tjfl 03 cC P . . 

p ^ -H XI >>rPi ® fiOH 
cs 3 a o tsD o > 
pH P • o 

H p ® -^ © H 
X O P> pop 
&E p 03 O p H 

-P += X 'H 

•H X 01 += p 
X O I 3x 

P CM a 03 O 

O - - . 

CM 


a 03 X 

HP 0 ) Pi 
05 P O X 03 rd 

P pii O X O 03 

b0& p 03 p p 
PQ>H-HOHpa 200 
(i3p03aooopox 
p p p p ph a 
a 3 o 5 PrCiHP 2 i ‘rtP 
rt xo3xp apprPi 
03 oXir\a 3 H-t^ a 30 x 
>j p O p p -P X -H 

5)0 03 • ax® Pi-P Xl > 

P Pi PH X 03 tiO 


03 Ogc3Q3a03aOOHp PiH 03 
op <1 >POPOChXc30^ ©P 
P>ch* oo^Ha ap f»0po5 
03 O H p p a ® ® >5 p p, 

X'H'HHO-HP .XPPPHP-Hxa 
03 03Hro03XXMpxaP 

ppti0t>>>ooa}®^xx pHoo • 
QJPClHHHapXO X-Hpp d 
P ® -H p cd X tjD WH H ® ox 03 cdHX 
oxH®Si®o-H-H pace oop® 
CP X PiH ®PXaXCMo5!X)OLP»PP®H 


p p 
® o 
H N 


o 
tiD 
p ® 
O P 
N O 

!>>x 


o H ';3- H a 


® 

® a 
P -H 
® O O H H 
P H 03 
a p X 

^ S’ll 
, o ® t(D 

H O X X -H 
O > K OH 
© 


pq p 
fctO 
H O 
P o 


p 

C ]3 H 

a P 


m 

m o>-^ 

O rt 

rt*— 03 

o 

VD 

Ch 

H 


-P o 
‘H ’H 

mxi -p 

'H 
'W. IS 


tjQ 
•H 

r-i 

■p Jh 03 
sd 0) g 
0) ‘H 
O Ch 
^ 'H 

0 4^ 

o a 

^ ‘T3 *H 
0^‘H -P 

CQ 
0 
H 

cd o 


o 

0 


03 


a 0 

4^ 

^ 0 
0 ^ 
H 


0 

H *H 

0 -p 

O -P 
> 

0^ 


a 

£53 03 X X O 

a 03 -H ® 02 

- p tio © P 

® -H p O 


© 

Xl 


p 

o 

•H * 
cd H 

CQ 0 
•H -P 
rH cd 

cd ^ 
P cd 
SU 


H EH O 


0 

0 


Td 0 

0 o 
tio 


cd 

t :3 0 

^ **0! . 

0 cp ch cd o 

^4^ O U U 

O 4^ 0 “ ‘ 

^ B U 


-p 

m 

Q 


0 


s:i 

. . . O rQ 
0 0 -H Ph 
j> <44 -P O -P 


P^•H 0 cd 0 c3 Is o 
5 h H < 1 > 
o IS H 
0 O 


H ^ 

0 • p 

0 ^ ^ 
H 0 
Pi cd -P F 4 
pi p 


0 
p p 
0 o 

^ <44 


ra Td 


0 0 0 0 €h 

0 Pig 0 O p 0 
cd o p p o o 0 

-P 4^ o p 0 -P 

p 0 p 0 

^ - 0 

4D 4 :) nd 
o cd 4 ^ p 


Pi cd Cd-P 
4 ^ 

p 


.g 


p 
o 
0 
H 

F^ 0 p 
o 0 Cd ^ 

O Pi O 

0 o a P» 

0 P o cd 


O4^Pi;dcdl0Pioa 


0 


0 P 


p a 0 

p 

P 3 C:^ 

0 

cd P >5 

-p 

rP 4 ^ H 

p 

0 0 cd 

p 


[>- 

H 


0 pp 

P cd P o 

0 cd o 

0 

0 0 H 

0 ^ 

p qd 4^ 0 P 

o p d cd 
p cd 


o 

VO 

Ch 

H 



4 ' 


ri TT a> 

0) OHS 
O 6J ft 

m -H a o 
ffi -p 3 c! 
ft-H tiD M _ 
tjD O 0) ,C| 

O -H O o 
O t:j a> o 3 
H — - (D 




ra I to ps 

I O -P H 

cs cS ^ a pi g 

a 5)D I 'H ft PJ rH r-( 

P4CQgpJP>4J-H» CD ft 

O-P-HCD pi O t>rt CD 

ftpsPia^-jOtao cDfflfciD a> ^ 

caftO-H^ Cc5ci5 HCDfj g+5 

Cj-P-Prd HiSft |SOtQ-H 

c3 pi 03 t :3 CD m >>-P H ccj a ‘Pi 

PlOpJHOPiPl isOCD f-lO 

4= HVOQpJOCD P(^fL(<DCD 

PiOtc>,y^(T3-H;El CD CDOP>CD 

p, CD O ^ P> -P p; |>>rO -ri CD P( 

aSNHrt'HOos CDPi*HOT:i+=> 

®>ia3-H taft c!5cq> pi 

PlHa CD© -HO H-dPi o© 

•H CC3 -H k5 O ^ H © 5 O CQ O 

HPiM+=>03+=dO spSoo 

oJcOaSft api PJO ft'© 

Jh a'draf^PiO 03,0 ‘Pi 

©ft plO Xt ©Oftft 

-POpl©H-HPlm OP>a-H*‘ p! 

p3 •H,c!a3H 05 o©-pH-p© 

ftpjcdftpift® a©pics©o5a 

a o -p oastj o-pofH^ tcD 

0-HPft-H©-H© Pi©^p50C\J© 

o -p o o m ra 64 -d ft d ftd t> ft m 


H 

iH 


<D 

0 


^ 

^ U U 1 


S4 Ci> 0 Q 

0 0 a 


<d tJ N “P 

xS O 


^ >s 

^ iH O 

U 

O O H Ph 

o O H 

0 

o 03 a 

o cd 'cJ -P 

00 CD O 

CO o ^ ^ 

P 

H U O 

iH P cd cd 

P4 

p ^ 



o >> 



{>^ p p 

U 


-p o -P 

0H 


CQ -P CQ 

^ cd < 3 ?^ 


(0^0 

0 o ^ 


rO 'H xS 

0-P p 4 cd 


O p 

g 


CO wM 



o 

<M 


VO 

VO 


cr\ 

cr\ 


H 

H 




48 


CQ -P 

■P O CQ 

rt -H © 

ffl -H -P ^ 

01 CS O O (!) 01 
p! 54 g ;C5 -P ?s 

H P" 3p rt o 

OQ © OJ <H -H -H 

© > CJ P 

•H © O 03 'H C3 

H O © -H O {> 

p msi 

s4iHC)8 ©©<i>ra 

© © ‘H © •© C? 

^ P 54 -tJ O -H O 
Ht(0©OS4s4O-H 
Cl3 "H ’H P cO © © P 

O ^ P © ^ 'CS © 

P W © 

M P 5:3 © !>>P 'd Q> 
©f3 -Hsai©Ps 3 P 4 
hq©PO^|:S©o 


© 

a 

•H 5 ^ 

P P O © 

© >s\ <H S 3 

© a 5 h op 

H © © a -H o 

© ft 14 l> 3 P 

•H 01 tjDP Sh © © 

p -H O tiOP H 01 


01 

I © 

O O p H 'H 

OJ P ft © 54 

© S 3 a 4 © 

© pra-HPOPp 

OP -HP •HOC©'© 


rl O ^ O P rQ 


.H(0|4 o©©544©o 


H © © ft O © 

j>-p 00 4© 01 fit 
© 54 P S3 

moi Q 


P©©!4HPf4!>POS3©-H 


O • O ft 54 O O 
© © © P © O © 
4 H 53< © O -© ft 


P P -H 

• O'© tiO'CS -H • 

B O © p p © ra-^ 
N p ©^ P © 
54ffl-H©4Sl4©® 


© © © prop 

p > a ^ p 4 

P 4'PP ©TOOPP 
ft©P44mrapo 
p p p © © © © 


© prop csip©^s© 
{10»>P 4ffl-H©4Sl4©® 

©©opp -HP rao©4 P 

4 mmpo -H a 4 ^H 

©© ©©■©'©>©>5 © 


Ppoo|>HP4ftq 


pTim©<j}o© 4 ®*'^ ••©^ 4 P^ 
5 ft i>p 43 o©-H^ 

OHra»P-HO©©4raH -©4 
4©© ffl©t(D t)Df404Ca©p© 
a?©P4Ht> PP *4© ©H-© 

+3 a .H © ft-H ©©40 p.oop 

H tiflIsatiD >®©tiDa 4 PpP 

•HP -HO© tH.Hi>©Op©H-HO 

phch-pph©©® £)0© anpp ft'-' 42 


iH 

© 

P 4 4^ 
p ® © tj >3 
© N P p3 © 
A >» P O H 

©H ft ft 

© a •© © 

CM P O p -H 

i-l P O P ©I 

P 

4 

P «© ® 
o p ft 
M © ft 
•H © 


P © •© 

p © ©P © 
© © p P 

O £(D o © -H 

4 P P O © 
© -H © 4 4 
ft p 4 © P 
•H © ft P 
O © P P 
O 4 P CMw 
H P PCI^ 


r© © 

© 4 

P © 

•H ^ 

a © © <© 

4 © © p 

© ft o © © 

p © p ® 

o © ©H 

rCS © © 4 © p 

© ©H © H a 

4 O © P © © 

ft © a p a <4 

Sj P 

0^0 so ■'d- 
iH ftCM 00 H H 

© 

» © -O 4 
O p H © © 

O 4H©H©P CQ 

W P©|>H4P 14 

© 0|>©00© P© 

4^P®^ 

©opft>s ©o pap 

!> © ftCQ i)Q >>H © © O tiD © 
©©+3 4P't31>|>‘H©ft 

© .©.HP -H P CO QQ 
© 40 - 0 © ©©'PPP. . 


■H P CO QQ 


© © P 43 O © 


©+344_© P44po 
>.«h© 0;PP'P,P©04 
4 © P P P P P '50-4 4 t)D 

© H © -H -H P 'H -H P rd p 

p © -H !s • © -H © -H 

pP4ra 4!>s:sP rop 
•HP -H P 4 H © H P 'H 
O -H iHI iH Cd CD 

OPH©H'c! P ©pP-dP 
©d pft dp-HOpen 
©OH® p©© ft p^ — ' 
a P ©•H'H|>'da 

4 |> t>j-H {> © © © -d P 
O © 4 tiO,P H © t>P © O H 
C\IP©4ft©'d -H 4 -HO 
HP®P'd© ©4'dOp}> 
P'HpP40©ra©5:3PPP 
<4 ch.H © fciOOdtO'd ffi ©P-H 


PH©C) 0 t> 
P ^ 4 P , 

O P ’H >3 

© H -d 4 


i 

H '-'N 

•H a 4 H a >s © 

«H © O ©M H 4 

P O ^ P © 

t!0©\ ©tr\ O P 
p >3pq H • H 

•H © "d tA ® -H 

ft .P"^ *504 'd 

P 4 JiDO p 0 

^ © OJ -H ITv P © 
KJP'-'j:3CM 4H P 

'd 

P P 

©©^ 

© © m 

©pft 
ft CQ'~^ 


C\J 

Os 


ved ) . 



49 


APPENDIX B 


MAIN 

PROGRAM 


SUBROUTINE 

LCPAS 


FUNCTION 




SUBROUTINE 

PRCIP 


jFUNCTICN 1 
! XMAX • 


^SUBROUTINE 

; POZB 


FUNCTION 

DIGFII 


SUBROUTINE! 
QUAD, J 


COMPUTER PROGRAM ORGANIZATION 


FUNCTION DIGFII - 
FUNCTION DIJFIL - 
SUBROUTINE LOPAS- 

SUBROUTINE PROLP - 
SUBROUTINE POZE - 


SUBROUTINE QUAD - 

FUNCTION XMAX - 


computes the low pass filtered sequence, 
computes the hand pass filtered sequence 
performs low pass filtering (for 500 Hz)on 
the speech input 

computes the coefficients of low pass filter 
computes the poles and zeros of the 
transfer function 

computes the quadratic expressions for 
transfer function 

computes the maximum value of an array. 









PLCI CHAET i.;F GOME^UTIR PROGRM 














THIS Pi\C<AGt FIMDS THK i NT ONA T ItTT PATTERN OF AM UTTERANCE 3F DUR/ T 
1. Z» SEE* 

u VECTOR SLI(l) contains SPEECH ir4PUI SAMPLES, 

DIMENSION ACS,E,3) ,E 16) ,SLI (2 S: ) , SGL I ( 5 E ) , SLO ( 3 ■ ) , VALR I 4 ) 
,VIN1(S,4) ,CNT(6) ,FR( 5) ,8(6,2, 5) 

INTEGER DfcCN(i0,40) 

COMMON ViN 
NTA3E=4 
PI = a.l.4.59 

c read speech sample, s i:rrc- slki) ■ rom tap,;:. 

REAOiNTAPE. )SLI 

C CALL LOPAS TO ELI Mi NaTt: FR - QU;E NC: ES ABOVi. HZ. 

CALL luPASCSLI) 

C DE-SIGM PROTOTYPE LOWPmSS FILT,,R 
u TR=TRaNSITI UM RaTIU. 

L ? B 1 = P S S B A I'i 0 i ■ j L . ' R A L c . 

t SBT^STDP oAiOD TOLcRAiTCE. 

L -SAK.=.SAMPLI NG FRCOUENCY. 

C' :UTFFL=CUT OFF FR.,,.QUcNCY OF PROTOTYPE LOWPAS FILTtRt THIS cQ, 

C BAND WIDTH OF THE C’P FILTER) I 

SATA TR, PBT,S8T,FSAM,CUTFFL/.9, .2, . 1 , 000. , 2 5. / I 

PISM=Pi/FSAF ' : 

w GLNcRATE FR;:,gUENCl t:S FOR 8P FILT R : 

FR{i)=7;. i 

D06SI = 2,... 

FR(I )=FR( 1-1)+CUTFFL 

o3 CONTINUE i: 

... SUBRJUriNE PROLP RETURNS CDEFFCIvNTS OF IP FILTER IN ARRAY A(I,J,Ki 

L I -INDEX OF Section cascade [ 

u J-IS : FOR NUMERATOR i 

L J-IS L FOR DENOMINATOR i 

C K-l,£,3 DESIGNATE COEFFI ECIENTS OF DELAYS | 

L. NSEC-NUMBLR OF SECTIONS IN CASCADE . I 

CALL PROLP{CUTFFL,TR,PBT,SeT,FSAM,E,A,NSEC ) S 

K = 1 [ 

u COMPUTATION BY FILTER-BANK : 

C NO. OF FILTeRS=EO ' 

D07JMFIL = l,i:: 

C FL-LOWER CUTOFF OF 8P FILTER 

u FH-JPPER CUTOFF OF BP FILTER [ 

FH=FR( NFIL+1) . j 

FL = FR(NFI.L) 

AL P= ( COS ( PISH* (FH+FL))) /COS tPiSM*(FH-FL) ) 

L COMPUTE COEFFICIENTS B{I,J,K) OF BP FILTER, AFTER APPLYING FREQUENC! 
t TRANSFORMATIONS TO LP FILTER. , | 

D025 1=1, NSEC - j 

D025J=l,2 ■ . 

Bd, J,i)=A( l,J,l ) : 

B( I, J,2)=ALP=«={A(I, J,2I-2.*A(I,J,I )) 


4 



B{ I, J, ; ) = {A (I, J,i )-A{ I, J,2) +A i I, J-,3) I , Jf 

B{ 1» J» 4}=ALP*!'{ hi I (I jJt» i ) ) 

B( It J»5}=A{I,J,3) 

5 CONTIMUE 

DOi5I = :.,ilS£C 
DO-ii^ J = .l t 2 
D035 <=2,5 

6(ItJtK)=B(IfJiK)/'’i(,jJ, ) 

cNTi n=B{i ,1,; )/B { i t , , i ) 

...5 CONTlNfUL 

D036I = i tS'jSEC 
DO 3 5 J= i 1 2 
B C I f J t , , ) = i • 

CONriNiUA 

c iNiTiA^is. viri(i,J) = :. 

0071 1 = ., , o 
D37iJ = i 
VIM« I, J )= .0 

71 CDMTIMUE 
MB = i 
N0=. 

DOTE IW’7aW=l,2-L 

C READ 57 SAMPLES OF INPUT INTO SG.. I ( I ) . 

0075 J= ■ t 5 ' 

SGLI ( J) = SLI (NB). 

NB=NJB+,, 

Vi CONTnUi: 

OUlk IRI;P = :., 

7 OUTPUT OF BP FILTER IS IN ARRA^Y SLO(I). 
i:)0?2tISAM=,. , 5 
X=SGLI (ISAM) 

SLa( ISAM)=DIJFIL(X,NSEC,B,E ) 

74 CONTIMUE 

C COMPUTc RMS VALUE OF INPUT AND OUTPUT SEyUHNCES. 
SUM1N=.. 

SUMOT=D.O 
D095ISQ=i, ,50 

SUMI N=SUMIN+S3LI ( ISQ)**Z 
SUM3T=SUM0T+SL0{ISa) **2 
95 CONTINUE 

C cVALUTE RATIO OF THE RMS VALUES OFTHE TWO SEQUENCES. 
RATIO=SQRT(SUMOT/SUMIN) 

VALRdWNDW )=RATIO 

C COMPARE RATIO WITH THRESHOLD VALUE THRS. 
IF(RATia.GT.THRS) 001094 
DECSl(NFIL,IWNDW+i)=. 

G0T0T2 

94 DSCN{SiFIL,IWNDW + l)=i 

72 CONTINUE 

PRINT 174, (VALR{IK),IK=I, 38) , 



io CQNriMU:- 

OUTPUT THL C-MATRiKj 

■ J-TlMt GCUMTERfQNE IWCREMEivT REPREGENTo 

^ 1- BP FILTER counter, 

PRi'^T Li? 

PR1MT2.. 9 
PRm?:;, 6 

D09ZIM = :. i,„U 
K = :,i-IM 
IL=K 

+ l ... , , 

PRl ^Ti f FR ( I L) 5 FR { i H ) » { DEC N{ R, J ) » J = I ? ) 

Vet CONTIMUi 

PRI'^Ti?."- 
PRIMT2?: 

,02 FOR'^ATC K, i.C.c,., 0.3 ) 

,,;,A FOR'1 Af ( A t „9F:.. ) 

OiS FQR'^AT t Xj ^FILTER B,:'iL'*/) 
wLO format t / ,&~Xs*TIME -*») 

, 03 F □ R '1 A T { i HU , 2 F 5 . ' » » 2 a i : ' 9 1 ) 

...9 FUR'^ATC (/)) , ^ ^ 

/2- FORMi^T t // /f 25K » * R «v C-M A T R I X FOR 

2 P A T T i_ R 

, ^,,5 FOR'^AKX, „9F&.^,tR’X,cF6.4 ) 

2,J FORMAT (i HI) 

STOP 

HMD 


53 MSEC. FOR KHZ. 


& N I iM T 3 N A T I 



«iSFTC 

fun: T I DM D IJFI L { X , Ni B , E ) 

C COMPUTES IHc BP FILTERcO OUTPUT OF TH£ INPUT SEQUENCE. 

u viNdfj)- Temporary storase for eequence 

L J- RfcPRGSEMTS DELAY NUMBER 

DIMENSION B(6, 2,5) ,b'{6) ,VIN(6s4) 

COMMON VIN 
ODlI=i,N 

P= { X-3 ( I , 2 , 2) 1 N { 1 , ) - C U , L, .. ) ^ Y I N{ I , ? ) -B ( I , £ , 4) N C I , 

i B { 1 , 2, 5 ) ^S'VI N{ I , N ) ) /B ( I , : , . ) 

X= ( 3 1 1 , 1 , 0 ) B ( i , , , ) =5= V I N { I , ) + B ( ! , 4- , 3 ) *V I M { 1 , 2 ) +B ( I , , 

2),+ B( , 3 )*Vi:H I ,4) ) *E( I )«2 .... 01692/0. 527 a97:, 

VINI i, x)=VIN(I ) 

VI N( i,E-)=VIN(I ,w) 

VIN( I,.:)=VIN(1 ) 

VIN{I,.)=P 

CONTINUE 

DlJFiL=X 

RETURN 

END 


E>- 

‘t)=xvIN 


I 



SUtH3JTINE PR3LP (0M:,,TR,R; ,R , ,3MS , E , A , ) 

COMPUTES THc COcFFiCICNTS OF THE CASCADED TRANSFER FUNCTiJM. 
DIMEMSiO.M „( 5),A(6 ,2 , 3 ) , G AMA { 21 ) , ZERO t 2.;. ) t UP ( A. ) , K P { ) 
CDM^LEK POLEC2 ) fGUM 
£P = SaRT{ (Ri*2. )~{R ^*2) )/(: .-MU 

AL = 5yRT(.,.-RZ*RU/(EP*RA ) 

0M2=Q^’ /TR 
XK = I ./(?:,4=OM5) 

C 3 M = T A N { 6 - 2 8 - . A * X K .1 .-i , , ) 

DON=i. /CON 
GON=D3N*DO!'; 

G A “ 4 L 
1 = 0 
1=U- i- 

GAMA {I )=GA 

G A = ( ^ ^ • ) / ( M ^ Q R ' { ) ) 

R 0 M i = { Ct M S / 3 . .1, ■:> 5 9 ) « T A .'4 { C D N* A A ) 

IF{RD^;..GT.0MU30TU . 

TR=3Mj /ROME 

NT = { )/ > 

CALL PDZ:UEP,GAHA,i» ZERO, UP, P) 

CALw aUA.D (DON, GulM,MT, ZEROtUP, '.P, c, A) 

J = 

XE = XMAXC ... ,NT) 

J = J 4 

lF(E(J).,rQ.X:. )GOTO : 

GOTO L 

EU J=E( JJ*{ l.-R. } 

RETJRM 

END 



^^'13 ft: 

SUB ^ OUT I POZ t: { E P I G AM A , NT t ZE P.0, U P , X P 3 

C COMPUTES THE POLES AND ZEROS OF THE TRANSFER FUNCTION. 
COMPLEX POLE {ZO), GUM 

DIMENSION GAMA(£l},Z.:R3{Ea) ,UP{EC),XP(2. ) 

X=].. /EP 

ZEROCi )=SQRT(( .+GAM (I) )/E.) 

POLE (i)=CMPLX( :.,X3 
X=ZER0{i) 

GUM=CMPLX{X, .) 

X=GAMA{i ) 

POLE a ) =GUM*CSURT( {CFPLT( . )+PDL£{ I ) ) / { CMPLX { X , 0 . J +P3L 
IF(N r..O.. )G;iT J 0 
DO 1 I=E,NT 
NJIG = { i:**l ) /T 
X = GAMA{ I ) 

GUM=CMPL.. (X, . ) 

DOI J = ; tXJIG 
NG=NJI G+J 

ZER3(NG)=SQRT{ ( I ' . + a)/ 2. ) *{ i. - ZERO { J ) )/(X-ZERO( J) ) ) 

ZERO! J > = SGRT(( { . + : .+ Zi.RO( J) ) /CX + ZERO{ J ) )) 

POLE ( MG) = CSGRTi ( CiMPL . ( t . , . . ) + GUM 3 =«=( CMPLX -POL£{ J ) ) 

, /(CMPLX( . . ,G. ) 'i'lGUM-POLEi J) ))) 

POLE! J )=C5QRT( ( CMPLXU. ,D. ) +GUM ) * ( CHPLX { i . , 0 . ) + POL fc ( J ) ) 

: / { C M Pu X { , f . ) * { GUM+ POL : ( J ) ) ) ) 

CONTINUE 
• N J 1 3 — { J., ^ N 1 ) / 2 
002 J = 

UP{J)=Rt AL{ POLc{ J) ) 

XP(J)=AIHAG(POLt ( J) ) 

CONTINUE 

RETURN 

END 



=i'iBFTC 

SUB^ OJ T I QUAD { DON, GOM , NT , ZERO, UP, XP , E , A } 

V COMPUTES THr QUADRATIC fcXPRESSIOri FOR FlLTcR TRANSFER FUNCi lO.M 
OIMENSIOM ZERDtZO) ,UP(2D) ,XP( 0 ) , c ( 6 ) , A ( 6 , 2 , 3 ) 

do:;: i=.,,nt 

A{ I, L,i )=GON+Z£RO{ 15^*2 
A ( I, 1, C )=-'£.=«'( GJN-2;,RU( I )*^ 2 ) 

YY=JP( 1 } ^=^2 + XP ( i ) 
zz^z.i'ttbsi xPd ) j 
A( 1,2,3 ) =YY-ZZ+GON 
A( I, 2, 2)=2.*{YY-G0; 0 
A( 1,2,. )=YY+ZZ + 30G 
S U ) = { A ( i , 3 , » } + A { I f - J 3 ) + A ( I ,2 . 3 ) ) 
i / { A I i , , ) +A {is: , • ; 3 + • { I , , , ) 3 

CDNTINJl. 

RcTJR'3 

END 



FUNCTi:-lNi XMAX{COL,!iJ 
DIM-: NS ION COL{6 5 
XMAK=3. 

DOil 1=1. ,N 

IF (( A3SI XMAX)-AbS (C3L( I ) )■) . GE. : , ) GO TO 
XMAX = 4t3S (COLII ) ) 

CONTINUE 

RETJRN 

END 



*IBFTC 

SUB^QJTINt LOPSiS(SLiJ 

SUBROUmt COHPUTcS TH£ L3WPASS i-ILTcRtD OUTPUT 
DIM£!^5IO''l SLK/: .T'0,C(5,2,3), C6),RINf6,2) 
COMMON} RIN 

DATA rRL,P8TL,S3TL,FSAML,CUTFFL/.8,.25» .i, 5,2 
DATA ./ 

CALL ?RaLP{CUTFFL,TRL,PBTL, SBl L, FSAML,E,C, N) 

ooi i=i ,£ :-j;: 

X=SL I{ I ) 

SLI{ I) =DIGFIL{ :\,N, :,C) 

CONTINUE 

RETURN} 

END 


♦ 


iij \J xu^ « f ‘ij m / 



-iBFTC 

fun:iidm OIGFIK V,N,„,4 ) 

DIM=iMSIO?4 a(&, 2,3) si:t6),RIN(6,2) 

COMMOM RIH 
DO i=I,N 

P={\//4 { I,c,i) )-{ aU , ,2 )/A( I, f .. ) )«RIN( I, 
( A{ 1 , 1 , 3 ) /A ( I , , 1 ) ) ^Ri Nt i , 2 ) 

V=£(l )*{A{I,., si 3*P + A{i,' ,2}=i'RIN(i,l) 

1 + A{I ,I i» 2 ) )^:,..19;>.692/. .3273-97. 0 

RIN< I, u)=RIN(I ) 

RIN{ I, i ) = P 

COMTI^UE 

DISFit=V 

REURvl 

EMD 



C THIS PROGRAM READS THE TAPE 
C TAPE T3 BE MOUNTED ON 8-4 
C ARRAY ADOUT CONTAINS ONE RECORD LENGTH 
C ARRAY DAT STORES EACH OF 12 BITS. 

C ARRAY PAT CONTAINS DECIMAL EQUIVALENT OF THEi2 BITS CORRESPONDI.N 
C ONE SAMPLE 

C IRGORD=NO. OF RECORDS TO B.-. READ 

DIMENSION PAT {4;.;/ ) 

INTEGER ADOUT{i333)»DAT{4 :u)yA,B,C,D,E 
IRCORD = i-: 

D020IM=l, IRCORO 

0011=1,1333 

ADOUT(n=C.O 

C READ THE TAPE 

CALL READ (ADOUT} 

IFLAG=: 

J= j. 

DQIO I =1,1335 

C SEPERATE THE 36 BITS INTO ...2 BITS EACH i 

A=ADOUT(I) 

IF{A.3E.O) G0TO5 ■ 

IFLAG=i. I 

A=-A : 

5 C=A/E : 

8 = A-C#E 
D=C/c 

C = C-D*t:; ^ 

IFdFLAG .NE. .) D=D+£/2 i 

DAT{J)=D i 

DAT(J+i)=C i 

DAT(J+2)=B I 

J=J+3 I 

IFLAG=0 f 

.0 CONTINUE I 

REWIND 1 i 

READ (I J OAT ! 

C CONVERSION TO DECIMAL i: 

DO i5K = X,400.- - f 

PAT( KJ =DAT{ K) 

PAT(K)=lU.=4'(PAT{K)/4 96.) ^ 

.5 CONTINUE ! 

PRINT 103, PAT 
Yu CONTINUE 

.00 FORMAT (SUX, IT/) ■ ! 

101 F0RMAT<5X,1C012) : 

.U3 FORMAT) 5X,20F6- 2) : 

STOP 
END 



*IBMAP COMV 

C MAP SUBROUTINE TO READ THE TAPE 



ENTRY 

READ 

READ 

SAVE 

i f 2 » A- 


CLA 

3,4 


STA 

10 


ENB 

=0 


TCOB 

* 


RDS 

T84 


RCHB 

10 


TRCB 

ERR 


ENS 

=0C'.;D.. 


RETURN 

READ 

10 

lORD 

**, ,i3 

TB4 

BOOL 

2224 

lRR 

TRA 

S. JXIT 


END 
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